第 29 卷 第 4 期 Eom 学 R Vol. 29, No. 4 
1986 年 11 月 ACTA ENTOMOLOGICA SINICA Nov., 1986 








Se OF RB EF dH RE EB ie Bj 
a FF EOS ih Ge Ot St 


(北京 农业 大 学 植保 系 ) 


J. Deltour P. Dagnelie 
《比利时 Gembloux 农学 院 ) 


摘要 ”北京 地 区 秋 白 荣 上 蚜虫 种 群 在 其 一 定 发 展 阶 段 的 空间 烙 局 是 可 以 用 Pearson UT 型 分 市 拟 合 的 。 
用 蒙特 卡 效 试验 研究 这 种 分 布 的 三 个 参数 一 一 平均 数 Re, 变异 系数 Ce Mia RA Cs 的 估计 误 沙 同样 本 
RE KR Ao e On ORC ATT DEPTH PR I ME PUR. TE ae, PL PRE A 
Bk mrt BC A EIE RAR MERRER ACT 596 的 情况 下 ,样本 容 证 为 50 fer HIRE nT IE GPS 


PWIA By 和 Ce ATUN Pir Ct TAL RA BA A 50 f. 
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一 、 做 为 菜 好 种 群 空 间 烙 局 模型 的 Pearson Ill 型 分 布 


rud op SU or RR, MERRE EI RIZR RAS REE, ORE 
数量 特征 不 稳定 ,变化 较 快 ;加 在 发 生 后 不 久 便 可 达到 相当 大 的 种 群 数量 ; 图样 本 中 个 体 
数 变 量 的 离散 程度 较 高 。 这 些 特点 是 在 兼 性 多 型 、 世 代 重 营 和 有 有 逆 好 迁 飞 等 生物 学 特点 
的 基础 上 形成 的 。 因 此 ,用 传统 的 Poisson 分 布 、 负 二 项 分 布 、Neyman A 型 分 布 ,Poisson- 
二 项 分 布 等 离散 型 概率 分 布 来 拟 合 沫 蚜 种 群 空间 抽样 数据 显然 是 行 不 通 的 。 棚 据 氏 endalt 
和 Stuart (1977) 提出 的 拟 合 复杂 观测 频率 分 布 的 三 种 主要 途径 , 沈 佐 锐 等 (1985) 初步 
证 明 , 利 用 Pearson 曲线 系统 可 以 描述 北京 地 区 秋 和 白菜 上 桃 蚜 Myzus persicae (Sulzer) FA 
BE WF Lipaphis erysimi (Kaltenbach) 混 生 种 群 的 空间 分 布 型 动态 ,而 其 中 II 型 曲线 则 
可 描述 菜 姬 种群 在 其 一 定 发 展 时 期 所 星 现 的 空间 分 布 状态 。 这 里 系 以 一 株 哑 为 一 个 样本 
单位 计数 每 株 菜 上 蚜虫 头 数 。 

* Pearson. HI 型 分 布 的 概率 密度 函数 为 
Jæ) 一 T (x — Ey eto (bez e <0) 


式 中 ay 84 5 为 原始 参数 ; æ, PB > 05 Mb=0, 成 为 T-5 fio l 
Pearson TIL 型 分 布 的 曲线 形状 又 可 由 三 个 数字 特征 决定 ， 数学 期望 Ex. PRB 


本 文 于 1985 年 4 月 收 到 。 
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Co 和 偏 态 系数 Cs; 故 三 者 亦 称 形状 参数 。 原 始 参 数 和 形状 参数 的 关系 是 





we oes — 2? 7. 5 page: (1 — 265 
d C3 d Ex» Cv» Cs mo BE (1 e) 
当 从 Pearson M Blo) 45 ik RER, RAEE E Ex, Co 和 Cs 的 估计 量 : 

Bx = m $i = x 

2XK; — 1». 
Pus [remm N—i 

by = Ien 

(N —3)- C» 


二 、 和 蒙特 卡 洛 试验 的 设计 思想 


有 害 生 物 单 种 种 群 的 研究 存在 着 两 个 方向 ,一 个 是 种 群 数 量 动态 ,一 个 是 种 群 空间 格 
局 。 就 前 者 而 言 ,人 们 关心 的 主要 是 对 Ex 的 估计 ,种 群 数量 或 密度 直接 关系 到 有 害 生物 
所 造成 的 经 济 损失 和 人 类 的 防治 策略 ;就 后 者 而 言 ， 为 了 研究 一 个 样本 是 否 可 用 Pearson 
II 型 分 布 拟 合 ,还 要 同时 考虑 对 Cv 和 Cs 的 估计 ,用 解析 的 方法 来 解决 这 个 问题 是 相 
当 困 难 的 。 然 而 这 恰 是 蒙特 卡 洛 方法 能 够 显示 其 优越 狂 的 场合 ， 即 使 对 Bx 的 估计 也 许 
存在 着 数据 变换 然后 用 公式 求 筑 的 方法 ， 但 那样 做 比较 适合 于 只 《有 一 个 或 两 个 参数 的 分 
布 ( 丁 岩 钦 ,1980; Ruesink，1980) ,而 对 Pearson M 型 分 布 来 说 , 仍 不 妨 用 蒙特 卡 洛 方法 
讨论 。 

蒙特 卡 洛 试验 的 原理 和 步 骏 

1. 首先 确定 试验 所 要 求 的 真 值 , 即 Pearson II 型 分 布 的 总 体 参数 Ex. Coo 和 Cso 
然后 在 计算 机 上 产生 Pearson I 型 分 布 变量 ， 这 相当 于 我 们 在 一 个 Pearson DI 型 分 布 
总 体 中 抽样 。 每 产生 一 个 变量 即 抽 得 一 个 样本 单位 ， 每 进行 一 次 试验 即 得 到 一 个 祥 本 。 
在 一 次 试验 中 所 产生 的 变量 总 数 则 相当 于 样本 容量 。 

2. 每 次 试验 结束 后 ， 计 算 每 个 总 体 参 数 的 样本 估计 量 。 在 参数 估计 的 三 种 常用 方 
法 一 一 矩 法 ,最 大 或 然 法 和 适 线 法 中 ,我 们 采用 矩 法 。. 

3. 重复 试验 ”所 重复 的 试验 次 数 应 使 样本 估计 量 的 平均 值 达到 试验 要 求 的 精度 。 

4. 设 计 不 同 的 样本 容量 水 平 “ 对 每 一 水 平 的 样本 容量 进行 了 上 述 三 步 后 ， 逐 个 计算 
三 个 参数 在 各 次 试验 中 的 估计 量 之 平均 值 ,并 与 真 值 比较 ,计算 相对 误差 。 

5. 以 相对 误差 评价 各 参数 估计 所 要 求 的 理论 抽样 数 。 

为 了 在 计算 机 上 产生 Pearson II 型 分 布 变 二 ,我 们 参考 了 华东 水 利 学 院 (1980) 介 绍 
的 框图 ;这 是 1973 年 Whittaker 提出 的 舍 选 抽样 方法 ,具有 较 高 的 抽样 效率 。 

米 个 试验 的 程序 用 TRS-80 微型 机 的 BASIC Level I 语言 写成 ,该 程序 所 需 的 均匀 
$y #5281 ER WLIO MEK RND 产生 。 

LUSH AGU BRE EASA, EERE AAE: a 所 输入 
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的 真 值 是 以 田间 实测 数据 为 基础 的 。 即 通过 适 线 法 用 Pearson UI 型 分 布 拟 合 实测 数据 。 
如 此 得 到 的 Ex, Co 和 Cs 估计 值 作为 蒙特 卡 洛 试验 的 输入 真 值 ; b. 由 计算 机 产生 的 
Pearson III. 型 分 布 变量 本 是 正 实数 ,我 们 作 了 整数 化 处 理 。 





三 、 蒙 特 卡 洛 法 本身 的 糖度 问题 . 
Vk 为 被 模拟 的 Pearson III 型 分 布 总 体 参 数 的 样本 佑 计 曙 ,这 恩 随 机 变量 ，En, om 

为 其 数学 期 望 与 方差 , 工 是 抽取 的 样本 数量 。 我 们 用 的 样本 平均 数 
m= tS (1) 


来 估计 En, ML oo 时 ， 不论? HSM, (my 一 En) BEA DH NCO, En/ 
VL)。 于 是 有 











Pm — Bal uf Li (2) 
“sy Af 2a 
在 95.5% BE KEL, 
VL, (3) 
on 
假定 公 许 误差 为 ep HERK 
oed, (4) 
En 
则 由 式 (3) 有 
> PUY cup, 
Lz( ; ) | Con (5) 
式 中 
Con En (6) 


为 人 的 变异 系数 (从 树 锋 等 ,1980)。 

在 开始 蒙特 卡 洛 试验 时 ， 我 们 并 不 知道 Con 的 值 ， 可 通过 下 面 两 个 阶段 性 试验 E 
决 : 

(1) 预备 阶段 试验 

设 工 一 50， 并 设 样本 容量 有 EN = 50, 100, 200, 300,400 五 个 水 平 。 给 定 总 体 
参数 后 , 即 进行 Pearson II 型 分 布 抽样 。 得 到 的 每 个 祥 本 用 和 矩 法 估计 Ex. Co 和 Co 这 
样 ,每 个 样本 容量 水 平 的 50 个 样本 都 可 给 出 这 三 个 估计 是 的 平均 数 和 变 虹 系数， 结果 列 
TOR lo 

根据 预备 阶段 试验 的 结果 ,可 粗略 估计 工 们 ， 即 先 以 Ve. Ve 和 Vs 28 Com 分 别 
代入 式 (5)， 并 令 。 一 5 求 出 相应 的 工 值 后 语 逐 当 放 大 , 列 人 入 20 

(2) 检验 阶段 试验 
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总 体 参 数 
50 
Ex, = 19.6692 ius 
Cvo = 0.8008 200 
Csa = 1.9263 300 
400 
50 

TESI 2 
Ex, = 596.221 tá 
Cv, = 0.7313 200 
Cs, — 1.8181 300 
400 





— 


* Segno 

















Cv Vu Cs Vs 
19.4926 | 0.0974 | 0.7785 | 0.1446 | 1.6353 | 0.4153 
19.3942 | 0.0738 | 0.7986 | 0.0903 | 1.8331 | 0.2604 
19.7592 | 0.0577 | 0.8077 | .0.0781 | 1.8724 | 0.2655 
19.5237 | 0.0519 |. 0.8065 .|:0.0565 | 1.8821 | 0.2204 
19.6273 | 0.0332 | 0.7992 | 0.0527 | 1.9055 | 0.2041 
593.424 | 0.0955 | 0.7267 | 0.1581 | 1.7146 | 0.4127 
601.352 | 0.0854 | 0.7278 | 0.1014 | 1.7825 | 0.3717 
597.416 | 0.0474 | 0.7299 | 0.0675 | 1.7120 | 0.2384 
598.542 0.0335 0.7420 0.0574 1.8628 0.2024 
596.127 | 0.0311 | 0.7317 | 0.0424 | 1.7333 | 0.1369 

















* dido Ex, Coy Cs 分 别 为 由 L= 50 个 样本 的 Ex, Co 和 C: 计算 得 到 的 平均 数 (mn), Vx, Vv 和 V, 
为 相应 的 Con 的 计算 值 。 


H2 5% 人 允许 误 状 下 的 工 值 














50 
100° 
200 
300 
400 
500 

















计算 值 估 定 值 计算 值 ee 计算 值 估 定 值 

50 15.2 50 33.45 70 275.95 300 
1000 8.71 40 13.05 50 108.49 250 
200 5.33 40 9.759 50 112.78 150 
300 4.31 30 5.108 40 77.722 100 
400 1.76 30 4.444 40 66.651 80 
50 14.6 50 39.99 80 272.51 300 
100 11.7 50 16.45 70 221.05 250 
: 200 3.59 30 - , 7.290 50 90.935 150 
300 1.80 30 5.272 40 65.545 100 
400 1.55 30 2.876 40 29.943 80 





250 
150 
100 
80 
70 











X3 “检验 阶段 试验 结果 之 一 例 (计算 机 打印 输出 ) 


Ex 


19,0447 
19.6878 
19.6416 
19.5050 
19.6204 
19.6061. 





SAMPLE 9.14 
(1) POPULATION PARAMETERS OF P-IIL DISTRIBUTION: 
Ex = 19.6692 Co = 0.8008 Cs = 1.9725 
(2) ORIGINAL AND SECONDARY PARAMETERS: 

AF = 1.0281 BI = 3.6985 BT. = 0.0644 NP =1 OP = 0.0281 


0.0827 
0.0524 
0.0409 
0.0419 





.. 0.0392: 











0.7938 
0.7976 
0.5026 
0.7972 
0.7997 






' 0.1330 





Vu 


0.1013 
0.0729 
0.0622 
0.0482 
0.0415 





* 对 于 Vx、Vv， 其 计算 出 的 < 们 显然 小 于 Ve 计算 的 。 位? 故 来 列 入 或 内 。 





Cs 


1.6291 
1.7451 
1.8779 
1.9390 
1.8970 
1.9310 








Vs 


0.3653 
0.2948 
0.2318 
0.2400 
0.1966 
0.1648 








e* 








4.2 
3.7 
3.8 
4.8 
4.4 
3.9 
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以 表 2 给 出 的 各 工 值 代替 预备 阶段 试验 中 的 L 一 50 ,重复 预备 阶段 试验 的 步骤 ， 然 
后 以 式 (5) 计 算 ce。 试验 至 少 重 复 两 次 。 胡 3 为 其 中 一 次 计算 机 输出 结果 。 

然后 我 们 就 检验 阶段 的 两 次 试验 结果 检验 所 取 的 世 值 是 否 合适 。 如 果 在 某 个 水 平 上 
某 参数 估计 的 相对 误差 大 于 5%, 则 这 个 水 平 的 LL 值 应 再 放大 一 些 。 结 果 表 明 , 我 们 这 里 
选取 的 值 都 能 满足 蒙特 卡 洛 试验 本 身 的 精度 要 求 。 





四 、 结 果 
各 样本 容量 水 平 的 参数 估计 相对 误差 可 用 下 式 计算 : 
IEn 一 Cont Lem LN x 100% 
式 中 ,1.96 J& a = 0.05 时 正 态 分 布 的 双 侧 分 位 数 ，mm 和 Sm 在 计算 时 分 别 以 检验 阶段 中 
两 次 试验 结果 合并 后 的 Ex, Co 和 Cs 及 Ve, Vo 和 Vs RER 

PELLE 和 Sym SEES 
其 中 Sm; 一 may + Cono 此 结果 列 于 表 40 


Ki 参数 估计 准确 度 与 样本 容量 的 关系 (9525 置 佑 水 平 ) 
e, ft or) 


er = max | 


























UN ie OE HT Ce SP Cs 
3,14 10,2 9,14 10.2 9,14 10,2 
50 1.42 1.43 4.17 2.48 20.41 16.00 
^- 100 1.03 1.20 2.09 2.09 13.70 > 11.24 
200 1.01 1.87 1.95 2.26 9.18 9.15 
300 0.92 2.02 1.74 2.27 7.50 7.62 
40D 0,93 1.1 1.40 2.48 7.08 9.83 
500 1.14 ”了 .41 0.97 0.51 5.07 6.38 





表 4 指 出 ,在 允许 误差 不 大 于 5% RIRN P ERRERA 50 IF, FA EY E BP 
地 估计 Ex 和 Co 值 了 ;但 对 于 估计 Cs 值 , 则 样本 容量 应 为 500 左右 。 


五 、 问 题 与 讨论 . 


BAR 4 SEL; 

1. 仅 就 佑 计 Ex 而 论 ,样本 容 景 还 可 更 小 些 。 

2. Pearson TH 型 分 布 总 体 人 参数 的 不 同 , 对 参数 估计 的 准确 度 亦 有 影响 。 

以 上 两 个 问题 应 在 进一步 的 试验 中 综合 研究 。 

3. 在 样本 容量 小 于 500 的 情况 下 , BRR TT IEA Cs 值 桐 对 误 薄 大 于 5%, ME 
更 为 精确 的 方法 , 如 适 线 法 ,但 这 至 少 要 拥有 一 台 徽 型 机 ( 沈 佐 锐 等 ,1985)。 

4. 从 树 锋 等 (1980) 在 用 这 特 卡 洛 试验 研究 了 Pearson M 型 分 布 釜 数 的 估计 方法 后 
指出 , 撼 法 \ 极 大 或 然 法 和 适 线 法 对 Ex, Co 和 Cs 的 佑 计 精 度 不 同 。 所 以 ， 当 我 们 为 
了 研究 用 适 乒 法 估计 仿 数 的 问题 而 进行 宜 特 卡 洛 试验 时 ， 上 面 介绍 中 用 到 短 法 的 地 . 辣 部 
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应 改 为 适 线 法 。 

5. 试验 速度 评价 以 表 3 为 例 , 进 行 完 表 中 全 部 计算 (从 EN 一 50 至 EN 一 500)， 其 
过 程 需 要 连续 进行 将 近 40 小 时 ,这 是 就 TRS-80 机 而 言 。 由 于 占 机 时 间 较 长 ， 就 限制 了 
蒙特 卡 洛 试验 中 一 些 最 优化 课题 的 研究 。 在 这 种 情况 下 ,最 好 应 用 大 ,中 型 计算 机 。o 
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STUDY ON SAMPLING FROM SIMULATED APHID POPULATION 
BY MONTE CARLO TEST 


Suen ZUO-RUI Guan ZHI-HE 


(Department of Plant Protection Beijing Agricultural University, Beijing) 


J. DeLrour P. DAGNELIE 


(Faculté des Sciences Agronomiques de l'Etat, 5800 Gembtoux, Belgique) 


For mixed populations of Myzus persicae (Sulzer) and Lipaphis erysimi (Kaltenbach) on 
Chinese cabbage in autumn in the area of Beijing, some of their spatial distributions can be si- 
mulated with the type Ul of Pearson curve system (P-IH distribution), whose variables can be 
generated in terms of the Monte Carlo method on TRS-80 microcomputer. The P-IIl distribu- 
tion is definited with three shape parameters: Ex, the mean; Cv, the coefficient of variation and 
Cs, the coefficient of skewness. The present paper deals with relation between’ acouracy of the 
parameter estimates and sample size from the P-III distribution population generated by the me- 
thod in Whittaker (1973). 

Let 1 denote the estimates of Ex, Cv and Cs; En, aq respectively mathematical expects and 
variances of the estimates; L the number of samples to enter into statistical calculation; and e 
the. pefmissible relative error, thus. 

Ul (228) Cen, 


where Cw = 23. 
En 
In order to ensure the precision of Monte Carto test, the L value should be determinated. So 
firstly, an exploratory test is made: 
Let L=50 and sample size has five levels, that is EN=50, 100, 200, 300, 400. At the end 
of the test, Cvy-values are estimated by the moment method and then L values are gained when 
e=5%. 
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Secondly, the L values are verified. 


Finally, relative errors in the parameter estimation, corresponding to the levels of sample 
size, are calculated by 


En — (mq £1.968m[A/ L 
= mes {En tune Bal Ly x 100% 


The result shows that in the case of e=5%, th: moment method can quite exactly estimate 
Ex and Cv even when EN=50; but for Cs, the method requires EN=500 at least. 








Key words Aphid population 
Spatial distribution 


Sampling———Monte Carlo Pearson curves 





